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CLUTTER-DOPPLER,  SPECTRAL  ANALYSIS 
FOR  A  SPACE-BASED  RADAR 


1.0  INTRODUCTION 


The  effect  of  strong  ionospheric  scintillation  on  the  total,  clutter-Doppler  spectrum  incurred 
by  a  monostatic  space-based  radar  (SBR)  system,  whose  orbit  is  above  the  ionosphere  and  whose 
carrier  frequency  is  between  100  and  1300  MHz,  is  analyzed.  This  work  is  a  continuation  of  the 
analysis  begun  in  (Mokole  and  Knepp  1991)  and  is  based  on  satellite  data  obtained  from  the  Defense 
Nuclear  Agency’s  Wideband  satellite  experiment  of  the  late  1970s  (Fremouw  et  al.  1978). 

In  addition  to  the  standard  moving-target-indicator  (MTI)  concept  of  designing  Doppler  filters 
to  detect  targets  of  unknown  velocity,  an  SBR’s  design  must  account  for  clutter  spread  caused  by 
ionospheric  scintillation,  platform  motion  relative  to  the  Earth’s  surface,  and  intrinsic  clutter  that 
consists  of  radar  instabilities  and  internal  motion  on  the  Earth’s  surface.  If  the  mean  Doppler  shift 
and  the  shape  of  the  clutter  spectrum  are  known,  Doppler  processors  can  maximize  the  detection 
of  moving  targets  (Andrews  1976;  Andrews  and  Gerlach  1989).  The  performance  of  optimal  filters 
is  heavily  dependent  on  the  width  of  the  clutter-Doppler  spectrum.  Therefore,  all  nonstochastic 
contributions  to  the  spectral  width  must  be  removed.  Consequently,  the  components  of  the  clutter- 
Doppler  spectrum  are  analyzed  to  determine  their  relative  significance.  One  goal  is  to  make  the 
minimum  detectable  velocity  (MDV)  as  small  as  possible  by  canceling  the  clutter.  To  get  the  precise 
relationship  between  the  spectral  spread  and  the  MDV,  one  must  consider  the  implementation  of 
the  clutter  canceller,  the  required  improvement  factor,  and  the  required  probabilities  of  detection 
and  false  alarm.  This  is  beyond  the  scope  of  this  report,  but  a  qualitative  indication  of  the  MDV 
may  be  obtained  by  comparing  the  relative  sizes  of  the  spectral  components.  Generally,  the  MDV 
increases  with  clutter  spread. 

The  clutter-Doppler  spectrum  is  the  result  of  three  contributions:  (1)  the  intrinsic  clutter 
that  is  backscattered  from  the  Earth’s  surface;  (2)  the  clutter  spectrum  resulting  from  ionospheric 
scintillation  experienced  by  the  signal  along  the  propagation  path;  and  (3)  the  spectrum  caused  by 
the  motion  of  the  satellite  relative  to  the  Earth  in  the  antenna  beam.  The  first  two  are  stochastic 
in  nature  and  are  assumed  to  be  statistically  independent,  while  the  third  contribution  is  strictly 
deterministic.  The  spectral  contributions  associated  with  the  random  processes  are  represented 
by  Gaussian  functions.  When  the  deterministic  contribution  is  also  represented  by  a  Gaussian 
function,  this  permits  the  use  of  a  more  tractable  formula  to  calculate  the  total  clutter  spread  for 
each  range  cell. 

Following  the  work  of  (Andrews  1976;  Knepp  and  Dana  1985),  a  formula  for  the  clutter-Doppler 
spectrum  is  derived  for  an  arbitrary  antenna  pattern.  This  formula  is  applied  for  a  Gaussian  pattern 
relative  to  certain  angular  coordinates  and  reduces  to  a  particularly  simple  form  for  a  range  cell, 
when  certain  conditions  apply.  Specifically,  the  functional  dependence  is  also  Gaussian, 


SclU)  = 


K 

—7= - exp 

y/liocL 


(/-/o)2V 

2 aCL  ) 


where 

+  a2PMi 


(1) 

(2) 
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K  is  a  constant,  /o  is  a  translation  in  frequency,  and  aj^,  age,  opm  are  the  spreads  caused  by 
intrinsic  backscatter,  scintillation  induced  by  the  ionosphere,  and  platform  motion,  respectively. 
Values  are  then  obtained  for  the  three  clutter  spreads  for  the  Wideband  satellite  pass  (ANCON 
08904)  that  experienced  the  most  severe  ionospheric  scintillation.  The  designation,  ANCON  08904, 
represents  the  data  set  that  was  received  at  the  Ancon  (Peru)  site  on  the  eighty-ninth  day  (089) 
of  1977  at  4  a.m.  (04).  The  results  are  then  compared  to  determine  the  relative  effects  of  each 
component.  The  spreads  attributable  to  scintillation  and  platform  motion  are  calculated  from 
analytical  formulae,  while  a  representative  value  of  the  intrinsic  clutter  is  specified  from  previous 
measurements.  Since  the  expression  for  age  is  easily  derived,  attention  is  focused  on  the  derivation 
of  a  formula  for  opm.  In  addition,  age  is  computed  for  all  Ancon  passes,  and  its  cumulative 
distribution  is  obtained. 

2.0  GEOMETRY  OF  COORDINATE  FRAMES 

For  this  problem,  the  position  of  the  satellite  in  space  is  referenced  relative  to  three  funda¬ 
mental,  orthogonal,  coordinate  systems  (6  Bate,  Ch.  2;  ADCOM  1977,  Section  2).  These  reference 
frames  differ  primarily  in  choice  of  origin  and  fundamental  reference  plane.  The  origins  are  the 
points  on  the  Earth’s  surface  (a  topocentric  system),  at  the  center  of  the  Earth  (a  geocentric 
system),  and  in  the  satellite  (a  vehicle-centered  system). 

The  geocentric  system,  usually  denoted  “Earth-Centered  Inertial,”  is  nonrotating.  Its  funda¬ 
mental  plane  is  the  equator,  the  X-axis  points  in  the  direction  of  the  vernal  equinox,  the  Z-axis 
points  toward  the  North  Pole,  and  the  F-axis  is  chosen  so  that  the  system  is  right-handed.  The 
unit  vectors  along  the  X-,  Y-,  Z- axes  respectively  are  I,  J,  K. 

Since  the  observations  of  the  Wideband  satellite  were  made  relative  to  two  ground  stations  at 
Ancon  and  Kwajalein,  a  topocentric  horizon  system  is  necessary.  In  particular,  the  location  of  the 
radar  is  the  origin.  The  coordinate  axes  are  denoted  Xh,  Yh ,  and  Z*,  with  the  positive  Xh.-  and 
Vh-axes  pointing  south  and  east,  respectively,  in  the  plane  that  is  tangent  to  the  Earth’s  surface, 
while  the  Z^-axis  points  radially  outward  from  the  Earth’s  center  (up).  The  associated  unit  vectors 
are  S,  E,  and  Z-  The  radar  sites  measured  the  range  and  direction  to  the  satellite.  “The  range  is 
simply  the  magnitude  of  the  vector  [-R5].  The  direction  ...  is  determined  by  two  angles  which 
can  be  picked  off  the  gimbal  axes  on  which  the  radar  antenna  [was]  mounted  [(Bate  et  al.  1971, 
p.84)].”  The  azimuth  angle  Az  was  measured  clockwise  from  the  negative  X/,-axis  (north),  and  the 
elevation  angle  El  was  measured  from  the  projection  of  the  radar  line-of-sight  onto  the  X/,F/i-plane 
to  the  line-of-sight  (— R5). 

The  satellite-fixed  frame  has  axes  pointing  in  the  nadir  (positive  z)  and  satellite- velocity  (pos¬ 
itive  x)  directions.  The  remaining  axis  ( y )  completes  the  right-handed  system.  When  discussing 
the  antenna  pattern,  spherical  coordinates  of  this  reference  frame  are  handy.  Define  the  spherical 
coordinates  ( r,<f>,9 )  by 


x  —  t  cos  9  sin  <f> 

y  =  r  sin  flsin  <j>  (3) 

z  =  r  cos  <f>. 


where  9  €  [0, 27r]  and  4>  €  [0, 7r] . 
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Fig  1(a)  —  The  geocentric,  topocentric, 
and  vehide-centered  coordinate  frames, 
ffi.  =  we  K  is  the  angular  velocity  of  the 
rotating  Earth.  This  figure  is  not  drawn 
to  scale. 


Zh(up)  Satellite 


Fig.  1(b)  —  The  topocentric 
coordinate  frame 


Fig.  1(c)  —  The  Earth-centered, 
inertial  (geocentric),  coordinate 
frame 
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3.0  CLUTTER  SPECTRUM 

An  MTI  system  must  remove  the  nonstochastic  contributions  (platform  motion)  to  the  clutter 
spread.  Two  systems  that  compensate  for  platform  motion  are  generally  used:  the  Time  Averaged 
Clutter  Coherent  Airborne  Radar  (TACCAR)  and  the  Displaced  Phase  Center  Antenna  (DPCA) 
systems.  TACCAR  is  a  phase-locked  system  that  was  designed  by  MIT  Lincoln  Laboratory  to 
remove  the  average  Doppler  shift  of  the  clutter,  whereas  the  DPCA  concept  uses  antenna  and  MTI 
techniques  to  remove  the  spectral  spreading  induced  by  platform  motion  (Andrews  1976;  Andrews 
1978). 

Figure  2  illustrates  the  effect  of  platform  motion  and  ionospheric  scintillation  on  the  clutter 
spectrum  that  is  observed  in  their  absence.  In  particular,  Fig  2(a)  depicts  a  relatively  narrow,  zero- 
mean  spectrum  caused  by  intrinsic  clutter  for  a  stationary  platform,  no  ionospheric  irregularities, 
and  no  ionospheric  motion.  The  combined  intrinsic-clutter  and  ionospheric  effects  are  visible  in 
Fig.  2(b),  where  the  scintillation-producing  irregularities  result  in  spectral  spreading.  Ionospheric 
motion  is  assumed  negligible;  thus  translation  of  the  spectral  peak  is  insignificant.  Figure  2(c) 
shows  the  translation  }d  and  additional  spreading  of  the  spectrum  from  platform  motion  and  the 
antenna  pattern. 


(a) 


0 


Fig.  2  —  (a)  Intrinsic  clutter 
spectrum  for  a  fixed 
platform  with  no  ionosphere, 

(b)  clutter  spectrum  that 
includes  the  effect  of  the 
ionosphere  for  a  stationary 
platform,  and 

(c)  clutter  spectrum  includes 
ionospheric  and  platform 
motion  effects 


First,  the  translational  effect  is  removed  by  applying  a  clutter-locking  technique  (Shrader  and 
Gregers- Hansen  1990;  Andrews  and  Gerlach  1989)  like  TACCAR.  This  effectively  centers  the  notch 
of  the  response  of  the  clutter  cancellation  processor  at  the  average  value  Jo  of  the  clutter  spectrum, 
which  is  necessary  for  good  performance  of  processors  like  MTI  and  DPCA. 

By  assuming  a  TACCAR-like  correction,  the  resulting  clutter  spectrum  is  derived  for  any 
clutter  patch  within  the  antenna  beam.  Results  are  then  derived  for  the  antenna’s  mainlobe,  which 
are  then  simplified  under  certain  assumptions. 

3.1  TACCAR  Correction 

Let  a  scatterer  on  the  Earth’s  surface  be  located  at  Ps  =  {Rs,4>s,0s)  relative  to  the  satellite- 
based  coordinate  frame;  let  R5  be  the  position  vector  to  the  scatterer;  let  Pa  =  (Ra,<t>a,9a)  be 
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the  coordinates  where  the  antenna’s  pointing  axis  (boresight)  intersects  the  Earth’s  surface;  let  h 
be  the  satellite’s  altitude;  let  a  be  the  grazing  angle;  and  let  vv  be  the  speed  of  the  satellite  (Fig 
3).  For  SBRs  above  the  ionosphere,  <f>a  G  [0,ir/2).  Since  the  Earth  is  rotating,  Ps  is  moving  with 
velocity  ve.  Hence  the  velocity  of  the  platform  relative  to  Ps  is  vr  =  vp  —  ve,  and  the  Doppler  shift 
of  the  scatterers  at  Ps  in  the  direction  R.s  is 


,  2 ft  (  \  2JtVr  2J*V* 

JD  =  -  Yr  •  R-5  =  -  cos  7r  ~  - COS  7s  = 

c  V  /  c  c 


2ftvr 


2/fUr 


2ftVr 


sin  <f>s  cos  $s  — 


2  ftVp 


sin  4>s  cos  6s,  (4) 


where  ft  is  the  transmission  frequency,  R-s  is  the  unit  vector  B s/Rsi  7s  is  the  angle  from  vp  to 
R5,  7r  is  the  angle  from  vr  to  Rs,  and  c  is  the  speed  of  propagation  (2.99792458  x  108  m/s).  For 
analytical  manipulations,  it  is  easier  to  use  7 s  instead  of  7 r.  This  is  done  in  the  remainder  of 
the  report.  However,  for  greater  accuracy,  the  equations  after  this  point  should  be  rederived  with 
cos7r  =  sin</>rcos#r,  where  (f>r  and  7r  are  the  spherical  angles  specified  by  Eq.(3)  for  the  vector 
vr.  The  goodness  of  the  approximations  in  Eq.(4)  at  the  altitudes  considered  here  is  discussed  in 
Section  4  and  Appendix  A. 


Fig.  3  —  The  geometry  of  a  scatterer  on  the  Earth’s  surface 
relative  to  the  vehide-centered  coordinate  frame 


If  A  is  the  ground  patch  illuminated  by  the  antenna’s  radiation  pattern  and  if  ft  is  the  corre¬ 
sponding  set  of  angles  the  average  Doppler  shift  over  the  ground  patch  is 

h  =  T^(A)ll  Mr'4,'<>)dA  (5) 

A 

A  common  practice  (Andrews  1976;  Knepp  and  Dana  1985)  is  to  integrate  over  the  projection 
ft  of  the  antenna  pattern  onto  a  reference  sphere  of  radius  Rq  about  the  origin  of  the  satellite- 
fixed  frame,  instead  of  integrating  over  the  surface  area  A.  This  inherently  assumes  that  fp  is 
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independent  of  Rs;  however,  fo  does  depend  indirectly  on  Rs  through  the  relative  speed  vr.  Since 
the  difference  between  the  minimum  and  maximum  ranges  to  the  illuminated  ground  patch  of  the 
mainbeam  can  be  as  large  as  540  km  (292  nm),  it  is  unclear  whether  this  approximation  is  accurate. 
To  determine  the  legitimacy  of  this  approximation,  one  must  derive  the  relationship  between  the 
coordinates  of  the  satellite-based  and  earth-centered-inertial  frames,  calculate  the  Jacobian  of  this 
transformation,  mathematically  describe  the  surface  patch,  and  evaluate  the  integrals  in  Eq.(5). 
Since  determining  the  validity  of  the  approximation  is  a  nontrivial  undertaking,  this  work  uses  it 
without  justification,  but  cautions  that  the  results  will  reflect  any  errors  inherent  in  its  usage. 

In  addition  to  the  preceding  assumptions,  suppose  that  the  clutter  over  small  solid  angles  is 
isotropic  and  homogeneous.  Therefore  (see  Appendix  B  for  details), 

=  xd^fi)  II  >°ia=l ^  (6) 

o 


where  d fl  =  sin  <f>  dcpdff.  The  TACCAR-like  correction  f'D  to  the  Doppler  shift  is  obtained  by 
subtracting  the  estimated  average  Doppler  from  the  Doppler  shift: 


2  ftv 

I'd  —  Id  ~  Id  =  — —  [cos  9S  sin  4>s  -  cos  9a  sin  <f>a\ , 


or  for  (<t>,9)  G  O, 


f'D(M 


2ftVr 

c 


[cos  0  sin  4>  -  cos  0a  sin  <t>a\ . 


(7) 


3.2  Analytical  Representation  of  Clutter  Spectrum 


The  received,  Doppler-corrected,  clutter  voltage  Vcl  for  a  monostatic  radar  is  analytically 
modeled  as 

VcL(t,r,4>,0)  =  V{r,(f>,0)AIN{t,<i),9)Asc(t,4>,0)exp(i2irtf'D((j),e)),  (8) 

where  t  is  the  local  clock  time  of  the  receiver  and  (r,  <j>,  9)  is  given  by  Eq.(3).  The  quantity  V 
represents  the  received  voltage  for  nonfluctuating  clutter  and  is  directly  proportional  to  G2/r2,  the 
ratio  of  the  two-way  antenna  gain  to  the  range  squared.  In  particular,  the  expression  V  =  CqG2 /r2 
is  used  in  the  remainder  of  the  report.  The  parameter  Co  is  the  constant  1  V-m2  so  that  the 
dimension  of  Co  It2  is  volts.  The  factors  A/^r  and  Ase  are  the  dimensionless  envelope  fluctuations 
caused  by  random  variations  of  the  surface  clutter  and  the  electron  density  along  the  two-way, 
ionospheric  propagation  path,  respectively.  Hence 

VcL(t,r,<(>,9)  =  C0~~~AiN(t,<j),8)Asc(t,(f>,e)exp(i2ntf'D(4>,9)),  (9) 

r£{<t>,9) 

The  random  processes  A/jv  and  Asc  result  from  entirely  different  physical  mechanisms  and,  conse¬ 
quently,  are  statistically  independent. 

The  autocorrelation  Rci  of  Vcl  >s 

Rcl(t,4>,0)  =  (VCL(t  +  T,4>,0)VSL(Ty<t>,9)).  (10) 

Thf  angular  brackets  denote  statistical  expectation,  and  the  superscripted  asterisk  denotes  complex 
conjugation.  Substituting  Eq.(9)  into  Eq.(10)  leads  to 

RcL(r,r,d>,0)  =  -^-jG(^,0)|4(A/yv(<  +  T)A}N(t)exp{2*iTf'D(<t>,0))Asc(t  +  r)A^(O).  (11) 
r  \<P,9) 
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Since  A /#  and  Asc  are  statistically  independent,  the  expectation  in  Eq.(ll)  becomes  a  product  of 
expectations.  Hence 


Rci(r,r,t,fi)  =  -^l^\G(<i>,e)\*(AIN(t  +  T)A}w(t)exP{2niTf,D(<j>,9)))(ASc(t  +  r)A^(t)) 

=  ■^^\G(4>J)\4^p(2TriTf'D(<f>,e))RIN(T)Rsc(T),  (12) 

where  iZ/w  and  Rsc  are  the  autocorrelations  of  the  intrinsic  and  propagation  envelopes,  respectively. 

The  Fourier  transform  of  Rcl  is  the  clutter  spectrum  Scl  the  direction  ( <j> ,  0)  of  the  antenna 

pattern: 


Scdf,<t>,0)  =  ^{^^|G(^fl)|4exp(2xtr/{)(d),0))fZ//v(r)/Z5C-(r)} 

=  ^^|G(^0)|4jr{exP(2WD(^^))^(r)}  *  t{Rsc(t)} 

=  G(<t>,0)\4SIN{f  -  f'D(<t>,0))  *  Ssc(f).  (13) 


The  asterisk  symbolizes  the  operation  of  convolution,  and  T  is  the  Fourier  transformation  operator. 
Although  an  asterisk  represents  two  mathematical  operations,  it  should  be  clear  from  the  context 
whether  an  asterisk  denotes  convolution  or  complex  conjugation.  Note  that  Scl  is  a  convolution 
of  the  clutter  spectrum  caused  by  scintillation  Ssc  and  the  translation  of  f'D  of  the  spectrum  from 
motion  on  the  Earth’s  surface  S hv. 

To  treat  the  total  clutter  spectrum  of  a  portion  of  the  antenna  beam,  in  particular,  the  surface 
subtended  by  the  solid  angle  corresponding  to  the  mainbeam  of  the  antenna  pattern,  the  transform 
representing  the  spectrum  for  a  given  direction  (Eq.(13))  must  be  integrated  over  the  appropriate 
angles.  The  spectrum  is  first  analyzed  for  incremental  subregions  of  A.  This  requires  nomenclature 
that  relates  an  incremental  region  to  the  satellite-fixed  reference  frame. 

The  mainbeam  ground  patch  (Fig.  4)  is  naturally  subdivided  by  range  cells  (radial)  corre¬ 
sponding  to  unambiguous  range  Ru  [=  c/(2/r)),  where  fT  is  the  pulse  repetition  frequency].  These 
unambiguous  range  cells  are  further  partitioned  into  Nr  subcells,  with 


Nr 


Ru 


Range  Resolution 


F2Ru1 

1 

—i 

.2/rT-cJ 

(14) 


The  double  brackets  denote  the  greatest  integer  function,  and  rc  is  the  compressed  pulsewidth. 

Since  the  range  does  not  change  appreciably  within  each  range  cell,  it  is  assumed  constant 
therein.  For  constant  range  r  and  arbitrary  0,  the  relationship  between  r  and  <f>  is 


r  =  Re 


h  -|-  Re 
Re 


cos  <j>  - 


(15) 


Only  angles  out  to  the  optical  horizon  are  considered  so  that  4>  lies  in  [0,  arcsin(/Ze/(/j  +  Re))j  •  The 
right  endpoint  is  the  angle  associated  with  the  optical  horizon  and  corresponds  to  a  zero  grazing 
angle  (a  =  0).  Equation  (15)  is  always  well-defined  for  this  interval. 
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Fig.  4  —  The  projection  of  an  antenna  pattern’s  mainbeam  onto 
the  Earth’s  surface.  The  solid  curve  within  the  ground  patch 
corresponds  to  the  antenna’s  boresight. 


The  number  Mo  of  unambiguous  range  intervals  covered  by  the  mainbeam’s  ground  patch 
depends  on  4>a  and  the  beamwidth  (A</>)  in  the  <f>  direction.  As  (j>a  and  A(f>  increase,  so  do  the 
number  of  unambiguous  intervals.  Suppose  that  M  and  M  +  Mo  —  1  are  natural  numbers  associated 
with  the  left  endpoints  of  the  first  and  last  unambiguous  range  intervals  of  the  ground  patch, 
respectively,  and  that  me  {l,...,Mo}  denotes  the  mth  of  these  unambiguous  intervals  (see  Fig. 
5(a)).  Therefore  [(M  +  m  -  1  )i?u,(M  +  m)Ru\  represents  the  mth  interval.  As  stated  earlier,  each 
interval  has  Nr  range  cells.  The  nth  range  cell  of  the  mth  unambiguous  interval  is  specified  by 

j^M  +  m  -  1  +  niV  -1)fiu,  ( 'm  +  m  -  1  +  . 

The  composition  of  the  mth  interval  is  depicted  in  Fig.  5(b).  Clearly  MoNr  range  cells  completely 
cover  the  mainbeam’s  ground  patch. 

Let  demarcate  the  antenna’s  mainbeam  beamwidth  in  the  9  direction  (azimuthal 

beamwidth).  Consider  the  clutter  spectrum  S ™  of  the  nth  range  cell  of  the  mth  unambiguous 
range  interval  associated  with  the  ground  patch.  To  obtain  S™,  the  transform  of  Eq.(13)  must  be 
integrated  over  the  surface 

At rm  =  |(r,<M)  :  r  G  [(m  +  m  -  1  +  —jy—)Ru,  (m  +  m-  1  +  ^")^]> 

and 

However,  on  each  range  cell,  r  (and  hence  <}>)  is  assumed  constant.  Arbitrarily  choose  the  constant 
value  of  r  to  be  the  midpoint  of  the  r-interval  of  Am,  and  denote  it  by  rm %■  The  corresponding 
value  4>mi  of  4>  can  be  obtained  from  Eq.(15)  by  solving  for  4>  and  substituting  for  r. 
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I - 1 -  l  - 1 - 1 -  l  - 1 - 1 

MRU  (M+l)Ru  (M+m-lJRy  (M+m)R„  (M+Mq-I)!^  (M+Mo)Ru 

(a) 

Ru 

Nr 

- - * - , 

I -  l  I  -  I -  l  - 1 

(M+m-l)Ru  (M+m-l+S^-)Ru  (M+m-l+^-JR,,  (M+m^ 

(b) 

Fig.  5  —  (a)  The  Mo  unambiguous  range  intervals  within  the  antenna’s  mainbeam,  and 
(b)  the  Nr  range  cells  of  the  m*  unambiguous  range  interval  in  the  mainbeam 


Under  these  assumptions,  the  integral  is  evaluated  over  the  set 
where  t(4> )  ~  r^.  Thus 

S%(f)  =  IJ  -  f'D(<f>,0))  *  Ssc(f)dil 

K 


Unit 


~  9i  r-  r 

rmn  J d> *  J 6‘ 


X  *  Sscif)^^  sin  ^wn  d<t>d,6 


-  4hn)  sin<^ 


f 

Je> 


\G{4>rrm,Q)\* 


X  SIN{f  -  foi^O))  *  Ssc(f)d9.  (16) 

Equation  (16)  expresses  S™  as  an  integral  of  a  convolution.  For  a  slightly  different  perspective, 
take  the  inverse  Fourier  transform  of  Eq.(16),  interchange  limits  of  integration,  and  apply  some 
transform  properties  to  get 

=  C°(^  ~  Rm(.T)RxW  f~  |G(A»,*)|4«tp(2x.> »))<»,  (17) 

'mn  J  8*^ 

which  is  the  autocorrelation  of  the  total  clutter  spectrum  incurred  by  the  surface  within  the  (m,n)th 
range  cell.  Clearly  it  consists  of  three  terms:  the  individual  autocorrelations  of  the  intrinsic  and 
propagation  envelopes  and  a  function  that  in  some  sense  adds  up  the  Doppler  effect  over  the  antenna 
beam.  Denote  the  product  of  the  integral  and  the  constants  in  Eq.(17)  by  Rpm(t).  Thus 

Co(<f>™  )  sin  <f>n 


prrtn 

nPM 


(r)  = 


f  IGitmmQy4  exp(2*iTf'd{<f>mn,0))  de,  (18) 
Je< 
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K%{r)  =  UE{T)K£{T)KRt{T).  (19) 

Upon  taking  the  Fourier  transform  of  the  preceding  equation  and  defining  Sp ^  to  be  the 
Fourier  transform  of  Rp\j,  one  arrives  at  the  general  expression 


S%{f)  =  (sft(f)  *  S&U))  * 


for  the  total  clutter  spectrum  of  the  (m,n)th  range  cell.  This  power  spectral  density  is  the  convo¬ 
lution  of  the  power  spectral  densities  of  the  backscatter  and  propagation  envelopes  and  a  function 
that  incorporates  the  deterministic  clutter  contribution  in  the  antenna’s  beam  from  the  Doppler 
associated  with  the  motion  of  the  satellite. 

Since  an  antenna  pattern  is  usually  referenced  relative  to  its  boresight,  introduce  two  new 
coordinates  4>  and  6  that  represent  the  position  relative  to  boresight  (<t>a,9a)-,  that  is, 

<t>  =  <f>  -  4>a  and  6  =  9  -  9a.  (21) 

In  Eq.(16),  make  the  change  of  variable  specified  by  Eq.(21)  to  obtain 

C'2(tbu  —  rbl  1  cin  rb 

crrax  /  _  '-'0  V'Crroi  V’mnJ  'rmn  I  I n.t  ±  aw\  errvn.  i  c  r?l  t  ±  n\\  .  omn  t  c\  jn  \ 


>e‘-ea 


4SK(f-F'D(4m,9))*s&(f)dO,  (22) 


where 


Fcfi&rmiQ)  ~  fbi't’mn  +  <t>a,  0  +  0a), 

G(4>mn,e)  =  G(4>rm  +  <i>a,9  +  oa).  (23) 

For  narrow  mainbeams,  the  translated  Doppler  correction  F'd  is  approximated  by  (Appendix  C) 

F'da(4>  rrai  7  6)  =  B4>mn  +  AO,  (2  4(a)) 

where 


A  =  -w\  s\n9as\n<j>a, 
B  -  w1  cos  6 a  cos  <f>a. 


(24(b)) 


For  an  SBR  with  a  low  pulse  repetition  frequency  (PRF),  only  the  clutter- Doppler  spectrum 
of  an  individual  range  cell  is  needed  in  developing  an  effective  signal  processor.  On  the  other 
hand,  a  high-PRF  SBR  simultaneously  receives  returns  from  multiple,  ambiguous  range  cells  in  the 
mainbeam.  Consequently,  one  must  consider  the  spectrum  generated  by  all  of  the  ambiguous  range 
cells  associated  with  a  given  range  over  the  entire  mainbeam  when  developing  the  signal  processor. 
In  this  report,  the  high-PRF  case  is  considered.  The  set  of  complex  clutter  returns  (voltages), 
which  are  obtained  by  integrating  over  each  range  cell,  are  assumed  to  be  statistically  independent 
and  to  have  zero  means.  Consequently,  the  corresponding  autocorrelations  R™  and  spectra  S™ 
are  additive  components  of  the  total  autocorrelation  and  total  spectrum  of  the  mainbeam;  that  is, 

Nr  Mo  Nr  Mo 

W  =  E{E™}  and  5ct(/)=E{E5^)}-  <25) 

n= 1  m  =  1  n=  1  m=l 

The  parenthetical  expressions  are  the  superpositions  of  the  autocorrelations  (RqL)  and  spectra 
(.%,),  respec'mvely,  of  the  nth  range  cells  of  each  unambiguous  range  interval  within  the  mainbeam. 
Equation  (25)  is  written  in  this  manner  because  a  receiver  cannot  distinguish  among  S}Al,  . . . ,  5^n; 
instead  it  lumps  them  together  as  for  each  n.  One  typically  says  that  the  clutter  from  each  of 
these  range  cells  folds  over  into  one  value,  as  though  the  summed  clutter  spectrum  comes  from  one 
range  cell.  The  contributions  from  the  sidelobes  are  ignored. 
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3.3  A  Special  Case 


Equation  (20)  has  a  rather  simple  form  and  facilitates  further  simplifications  for  some  functions. 
In  particular,  suppose  S™  is  Gaussian  with  zero-mean  and  standard  deviation  ap  for  f3  G  {IN,SC}\ 
that  is, 


*r(/)  =  ^=-e* p  -O  • 


p  w 7  y/2 nap 

Therefore,  the  convolution  of  S iff  and  S^c  becomes 


exp  - 


+  ct5c)  ) 


Hence 


smn  = 


2(az/v  +  a\c) 


<Usr)) 


*«(/)• 


Upon  taking  the  Fourier  transform  of  Eq.(  ,  with  the  change  of  variable  given  by  Eq.(21), 
observe  that 


S^f(f)=  r  r m"  C°2^-.^)sin  ^  I  Qi^et  exp(27T  nF'oi^e))  exp(-2;r  ifr)  dUr 

=  r  *6{f-F'D(*n«j))d9,  (28) 

rwi 

where  <5  is  the  Dirac  delta,  =  0^  -  0a  <  0,  and  =  9 -  8a  >  0.  Formally  substituting 
Eq.(28)  into  Eq.(27),  interchanging  the  order  of  integration,  and  applying  the  sifting  property  of 
the  Dirac  delta  yield 

S%U)  -  f~  |g{^.g)l4exp  (-  (/  T  )  S.  (29) 

rlmV2irs/(T2IN  +  eric  \  is  +  °sc)  / 

It  is  tacitly  assumed  that  the  mathematical  operations  resulting  in  Eq.(29)  are  valid  over  the  class 
of  functions  that  are  considered. 

An  explicit  representation  of  the  antenna  pattern  is  required  for  there  to  be  any  chance  of 
getting  a  closed  form  expression  for  S^£.  Even  when  the  pattern  is  specified,  a  closed-form  expres¬ 
sion  is  difficult  to  obtain.  Usually  one  must  resort  to  a  numerical  evaluation  of  Eq.(29).  However, 
for  a  Gaussian-shaped  antenna  pattern  having  circles  in  the  (<f>, ^-coordinate  frame,  an  explicit 
representation  of  S™  is  obtained;  specifically,  it  is  now  shown  that  Eq.(29)  reduces  to  Eq.(l),  when 
a  certain  approximation  is  valid. 

Let  the  one-way,  power  gain  of  the  antenna  be  characterized  by 


\G(<f>mn,0)\4  e\p(2TriTF,D((j)rrm,9))  exp(-27r ifr)  dOdr 


1  G(4>,  *)|2  =  go  exp  J :  (30) 

where  go  is  the  power  gain  at  boresight,  and  r?  is  a  constant  that  is  chosen  so  that  the  power  gain 
within  the  mainlobe  is  95%  of  the  total,  two-way  gain.  For  narrow  mainbeams,  Fp  is  approximated 
by  B<t>mn  -I-  AS  (Eq.(24)).  In  this  case,  the  level  curves  of  the  two-way  gain  \Q\A  are  circles  in  the 


11 


MOKOLE 


(4>,  0)-coordinate  frame.  Note  that  these  level  curves  are  not  the  cross  sections  of  the  antenna 
pattern  that  are  perpendicular  to  the  boresight.  Define  2po  to  be  the  desired,  mainlobe  beamwidth. 
Since  the  Gaussian  pattern  has  no  nulls,  interpret  this  to  mean  that  95%  of  the  power  is  in  the 
region  0  <  4>2  +  62  <  pi,  which  implies  that  77  =  po\/8/ln(20)  (see  Appendix  D). 

With  the  Gaussian  power  gain  and  a  narrow  mainbeam,  Eq.(29)  reduces  to 

SZV)  -  2 [erf(C.)  -  erf(0]  ^  exp  -  (31) 


where 


\A\rj  T]  2ftvr ,  .  /,  i  .  j 
oPM  =  — p  — —  |sin0a|sm<£a, 

2  2  2  i  2 

°CL  ~  alN  T  aSC  +  aPM » 

rx  .2  VV*  Cli*™  ~  <tLi)  sin  ^  t  \ 


V  = 

lmn 


4<tcl 


+ ’so) 


4<t cl 


'IN 


+  <Tsc) 


-1?“  - 
v  mn 


At)2  (/  - <£„„) 


2  tx 

erf(z)  =  — I  exp(— 1 2)dx. 

V*  Jo 

When  the  integral  in  Eq.(29)  closely  approximates 

r  ,ca  „„ (_  u-nd™A)\  dS 

expl  2(^+0  )M’ 


(32) 


(33) 


Sq£  is  approximately 


s&Af)  = 


Kr, 


y/2n<T  cl 


exp 


(- 


CL 


(34) 


Therefore,  when  this  approximation  is  valid,  the  clutter-Doppler  spectrum  within  a  range  cell  is 
characterized  by  an  expression  having  the  form  of  Eq.(l). 

Clearly,  the  approximation  is  not  good  for  range  cells  at  the  heel  and  toe  of  the  mainbeam’s 
footprint,  that  is,  for  4>rm  —  ±Po-  For  further  details,  see  Appendix  D.  Since  Eq.(34)  is  simpler 
than  Eq.(31),  it  is  more  desirable  for  analytical  purposes;  however,  it  may  result  in  errors,  and  one 
should  be  cognizant  of  this. 


4.0  CLUTTER  SPECTRUM  FROM  WIDEBAND  DATA 


The  1977  Wideband  data  (Fremouw  et  al.  1978;  Knepp  and  Mokole  1991)  are  used  to  extrap¬ 
olate  the  clutter-Doppler  spectrum  and  spreads  for  a  sea-clutter  patch  and  for  the  Gaussian  power 
gain  of  Eq.(31).  Although  these  data  correspond  to  severe,  ionospheric-scintillation  conditions  for 
1977,  the  clutter-Doppler  spread  caused  by  scintillation  is  not  the  worst  (largest),  since  the  data 
were  measured  at  the  end  of  a  period  of  solar  minimum. 
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Data  from  22  satellite  passes  are  used  to  construct  a  distribution  of  <jcl-  In  particular,  a  cl 
is  computed  for  the  satellite  pass  (ANCON  08904)  that  experienced  the  most  intense  scintillation; 
that  is,  the  values  of  tq  are  the  smallest  (Mokole  and  Knepp  1991;  Knepp  and  Mokole  1991).  The 
Ancon  receiver  site  was  located  at  -11.767°  latitude  and  -77.500°  longitude.  The  orbit  of  ANCON 
08904  is  chosen  as  representative  of  a  retrograde,  circular,  nearly  polar  orbit  at  an  altitude  of 
approximately  1030  km.  The  position  of  this  pass,  measured  from  Ancon,  allows  estimation  of  vT, 
4>a,  and  0a. 

The  total  clutter- Doppler  spread  depends  on  three  parameters:  gjn,  use,  <?pm-  A  typical 
value  for  worst-case,  clutter-Doppler  spread  a  in  for  internal  motion  of  the  sea  is  1  m/s  (Shrader 
and  Gregers-Hansen  1990,  Table  15.1,  p.15.9),  whereas  the  scintillation  spread  a  sc  is  determined 
from  an  analytical  formula  that  depends  on  tq.  According  to  Eq.(34),  <?pm  is  calculable  once 
the  antenna’s  beamwidth  2po  and  the  parameters  vr ,  <f>a,  0a  are  specified.  This  report  considers 
narrow  mainbeams  on  the  order  of  1°  to  2°;  so  a  nominal  value  of  po  =  1°  is  selected.  The 
individual  clutter-Doppler  spreads  are  compared,  and  their  importance  relative  to  ocl  is  discussed. 
In  addition,  examples  are  considered  to  determine  whether  Eq.(34),  the  approximation  of  Eq.(31), 
may  be  used  to  calculate  the  total,  clutter-Doppler  spectrum. 

4.1  Clutter-Doppler  Spread  Caused  by  Platform  Motion 

Since  a  pm  depends  on  <j>a,  0a,  vT,  and  since  these  parameters  depend  on  a  satellite’s  orbit, 
orbital  data  pertinent  to  ANCON  08904  is  discussed.  The  satellite  was  launched  by  the  P76-5 
Scout  Rocket  on  22  May  1976.  The  orbital  elements  and  associated  information  (Table  1)  are 
determinable  from  a  North  American  Aerospace  Defense  Command  (NORAD)  pass  planner  and 
(King-Hele  et  al.  1988).  The  orbital  elements  of  Table  1  are  used  to  create  an  input  orbit  for  the 
most  current  version  of  the  Fast  Orbital  Prediction  Program  (FOPP8)  (Eisele  and  Shannon  1975) 
that  most  closely  matches  the  NORAD  orbital  data.  With  this  facsimile  orbit,  FOPP8  calculates 
<f>a,  &a,  Ra •  Then  the  range  outputs  ( Ra )  of  FOPP8  are  used  to  estimate  vr,  which  is  used  to 
calculate  a  pm- 


Table  1  —  Orbital  Elements  of  Wideband  Satellite 


SSDC  No. 

08860 

International  No. 

1976-47A 

Launch  Date 

22.32  May  1976 

Expiration  Date 

400  yr 

Payload 

P76-5 

Inclination  Angle 

1.73953751649  rad 

Eccentricity 

0.054830209 

Angle  of  Right  Ascension 

6.27122910366  rad 

Angle  of  Perigee  from  Ascending  Node 

1.60201995295  rad 

Normalized  Semi-Major  Axis 

1.1608806 

Mean  Motion 

0.0594559435781 

Greenwich  Hour  @  Epoch 

2.57899235231  min 

Period 

105.5  min 

Perigee  Height 

996  km 

Apogee  Height 

1060  km 

The  time  of  each  data  point  (Table  2)  is  the  local  mean  solar  time  at  the  Greenwich  meridian, 
variously  called  Greenwich  Mean  Time  (GMT),  Universal  Time  (UT),  or  Zulu  (Z)  time  (Bate  et  al. 
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1971).  The  entry  4:4:30.0  means  4  hr,  4  min,  and  30  sec  in  the  morning  at  Greenwich.  The  azimuthal 
Az  and  elevation  El  angles  are  measured  in  the  topocentric  frame  at  Ancon  (Fig.  1(b)),  while  the 
latitude  (Lat)  and  longitude  (Lon)  are  measured  with  respect  to  the  rotating,  Earth- centered  frame 
with  0°-longitudinal  meridian  passing  through  Greenwich.  From  the  latitudes  of  Table  2,  a  north- 
to-south  equatorial  crossing  occurs  between  4:6:30.0  and  4:7:30.0.  To  make  use  of  the  data,  assume 
that  the  mainbeam  of  a  hypothetical  antenna  points  at  Ancon;  hence  Rs  =  Ra ■  The  ranges  Ra 
generated  by  FOPP8  are  not  exactly  those  of  the  NORAD  pass  planners,  however  they  are  close 
enough  to  permit  calculation  of  vT.  The  angles  <t>a  and  9a  correspond  to  the  parameters  Dip  Angle 
+  7t/2  and  BackAz,  respectively,  of  FOPP8. 

Table  2  —  Orbital  Data  Relative  To  Ancon 


UT 

(hr:min:s) 

Az 

(Deg) 

El 

(Deg) 

Lat(+N) 

(Deg) 

Lon(+E) 

(Deg) 

Ra 

(km) 

{ Pa 

(Deg) 

0a 

(Deg) 

4:04:30.0 

23.6 

10.4 

8.1 

-68.6 

2822.4480 

58.297 

346.34 

4:05:30.0 

25.5 

15.8 

4.7 

-69.4 

2453.3444 

56.503 

344.68 

4:06:30.0 

28.1 

22.4 

1.3 

-70.2 

2097.0196 

53.516 

342.25 

4:07:30.0 

32.3 

30.9 

-2.0 

-71.1 

1762.7336 

48.752 

338.48 

4:08:30.0 

39.8 

42.0 

-5.4 

-71.9 

1467.5248 

41.427 

331.68 

4:09:30.0 

56.1 

55.7 

-8.8 

-72.7 

1240.6548 

31.188 

317.52 

4:10:30.0 

96.6 

65.8 

-12.2 

-73.5 

1125.6456 

21.731 

282.61 

4:11:00.0 

124.1 

64.7 

-13.8 

-74.0 

1156.2036 

24.411 

232.50 

4:12:00.0 

157.9 

52.8 

-17.2 

-74.8 

1322.3280 

35.007 

207.69 

4:13:00.0 

171.3 

39.7 

-20.6 

-75.7 

1581.2376 

44.066 

197.24 

4:14:00.0 

177.8 

29.3 

-23.9 

-76.6 

1894.5960 

50.231 

191.86 

4:15:00.0 

-178.5 

21.4 

-27.2 

-77.6 

2238.6976 

54.182 

188.64 

4:16:00.0 

-176.2 

15.1 

-30.6 

-78.6 

2600.3932 

56.632 

186.50 

Since  Ra  is  the  distance  measured  in  the  topocentric  frame  at  Ancon,  Ra  is  the  relative  radial 
velocity  vr.  According  to  Eq.(A17)  of  Appendix  A, 

Yr  •  Rs  =  Yr  •  Ra  =  -Ra, 

where  Ra  is  the  amplitude  of  Ra.  However,  the  approximation 

vr  •  Ra  =  tv  cos  7r  ~  vr  cos  70  =  tv  sin  <pa  cos  6a 

implies  that 

tv  —  — - — — —  ,  (35) 

sin  4>a  cos  9a 

provided  the  denominator  is  not  zero.  Therefore,  to  estimate  the  speed  vT  at  a  particular  time  to 
of  the  orbit,  the  following  difference  scheme  is  employed, 

In/ONl  I  f  ■  ±  I  dl\l  Rai^O  A  — )  ~  Rai^O  2”)  /oc\ 

|Ra(<o)|  =  |vr  -Ra|  =  (*vsin<£a|cos0o|)|t=to  ~  - - —  ,  (36) 
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where  a  time  increment  At  of  2.4  s  is  chosen.  Substituting  At  and  Ra  from  Table  3  into  Eq.(36) 
yields  the  estimates  of  vr.  These  speeds  are  listed  in  Table  4  with  the  other  parameters  that  are 
necessary  in  computing  ap\i.  By  Eqs. (24(b)),  (32),  and  (35), 


OPM 


7T  /  8 

“  720  V  In~20 


8  2ft\Ra\ 


tan  9  n 


(37) 


if  cos 0a  ^  0.  The  evaluations  of  Eq.(37)  for  137.6748  MHz  (VHF),  413.0244  MHz  (UHF),  and 
1239.0730  MHz  (L-band)  axe  also  provided  by  Table  4.  Clearly,  crPM  increases  with  increasing 
frequency. 


Table  3  —  Ranges  and  Universal  Times  For  Ancon 


UT 

(hr:min:s) 

Ra 

(km) 

UT 

(hr:min:s) 

Ra 

(km) 

UT 

(hr:min:s) 

Ra 

(km) 

4:04:28.8 

2829.8560 

4:09:28.8 

1244.3588 

4:12:58.8 

1575.4964 

4:04:30.0 

2822.4480 

4:09:30.0 

1240.6548 

4:13:00.0 

1581.2376 

4:04:31.2 

2814.8548 

4:09:31.2 

1237.1360 

4:13:01.2 

1587.1640 

4:05:28.8 

2460.5672 

4:10:28.8 

1126.5716 

4:13:58.8 

1887.9288 

4:05:30.0 

2453.3444 

4:10:30.0 

1125.6456 

4:14:00.0 

1894.5960 

4:05:31.2 

2445.9364 

4:10:31.2 

1124.7196 

4:14:01.2 

1901.2632 

4:06:28.8 

2103.8720 

4:10:58.8 

1154.1664 

4:14:58.8 

2231.6600 

4:06:30.0 

2097.0196 

4:11:00.0 

1156.2036 

4:15:00.0 

2238.6976 

4:06:31.2 

2089.9820 

4:11:01.2 

1158.2408 

4:15:01.2 

2245.7352 

4:07:28.8 

1769.2156 

4:11:58.8 

1317.8832 

4:15:58.8 

2592.9852 

4:07:30.0 

1762.7336 

4:12:00.0 

1322.3280 

4:16:00.0 

2600.3932 

4:07:31.2 

1756.4368 

4:12:01.2 

1326.7728 

4:16:01.2 

2607.6160 

4:08:28.8 

1472.8956 

4:08:30.0 

1467.5248 

4:08:31.2 

1462.1540 

Table  4  Clutter- Doppler  Spread  From  Platform  Motion  at  Ancon 


4^<x 

(deg) 

9a 

(deg) 

1  Ra  | 
(km/s) 

vr 

(km/s) 

<rm@VHF 

(Hz) 

<rm@UHF 

(Hz) 

<TpM@L-band 

(Hz) 

58.297 

346.34 

6.2505 

7.5606 

9.95 

29.85 

89.54 

56.503 

344.68 

6.0962 

7.5797 

10.94 

32.81 

98.43 

53.516 

342.25 

5.7875 

7.5580 

12.13 

36.40 

109.19 

48.752 

338.48 

5.3245 

7.6124 

13.75 

41.25 

123.75 

41.427 

331.68 

4.4757 

7.6840 

15.80 

47.39 

142.16 

31.188 

317.52 

3.0095 

7.8799 

18.05 

54.14 

162.43 

21.731 

282.61 

0.7717 

9.5471 

22.59 

67.77 

203.30 

24.411 

232.50 

1.6977 

6.7479 

14.49 

43.47 

130.40 

35.007 

207.69 

3.7040 

7.2917 

12.73 

38.19 

114.57 

44.066 

197.24 

4.8615 

7.3189 

9.88 

29.64 

88.92 

50.231 

191.86 

5.5560 

7.3861 

7.64 

22.92 

68.77 

54.182 

188.64 

5.8647 

7.3155 

5.84 

17.51 

52.52 

56.632 

186.50 

6.0962 

7.3467 

4.55 

13.65 

40.94 
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4.2  Computation  of  asc 

When  severe  scintillation  conditions  exist,  the  autocorrelation  Rsc  (Eqs.  (11)  and  (12))  of  the 
envelope  fluctuation  Asc ,  which  are  caused  by  ionospheric  scintillation,  is  Gaussian  (Mokole  and 
Knepp  1991);  that  is, 


Rsc(r)  =  (. Aac{t  +  T,<t>,9)A'x{ty4>,9 ))  =  e"^. 

The  parameter  To  is  called  the  complex  decorrelation  time.  If  propagating  pulses  are  separated  in 
time  by  an  amount  less  than  tq ,  they  are  presumed  coherent.  It  is  assumed  that  Rsc  is  constant 
with  respect  to  small  changes  in  <f>  and  0.  Since  the  line  of  sight  between  the  Wideband  satellite 
and  Ancon  changes  slowly  over  the  orbit,  this  assumption  holds.  Equation  (37)  implies  that 

Sic(/)  =  iexp(^)'  (38) 

Clearly,  Ssc  spreads  (asc  increases)  as  ro  decreases.  For  severe  scintillation,  To  <  0.25  s. 
Figure  6  provides  histograms  of  the  measured  values  of  tq  for  all  the  Ancon  satellite  passes  at  VHF, 
UHF,  and  L-band.  As  the  frequency  increases,  the  percentage  of  To  at  the  smaller  values  decreases. 
In  particular,  slightly  more  than  23%  of  the  values  of  ro  are  below  0.25  s  at  VHF.  In  contrast,  less 
than  3.5%  and  0.11%  of  the  r0  are  below  0.25  s  at  UHF  and  L-band,  respectively.  Although  not 
shown  in  Fig.  6,  17.45%,  44.81%,  and  70.76%  of  the  ro  are  greater  than  8.392  s  at  VHF,  UHF, 
and  L-band.  Additionally,  the  minimum  value  for  each  frequency  is  0.024  s.  The  measured  minima 
are  not  exact  but  are  limited  by  the  integration  time  of  the  Wideband  receiver  and  the  sampling 
interval  (2  ms).  Hence  24  ms  is  probably  not  the  actual  minimum  r0;  however,  values  less  than  24 
ms  occur  infrequently,  since  this  value  is  not  measured  often. 

Figure  7  displays  the  distribution  of  (75c,  corresponding  to  the  measured  values  of  To,  for 
each  frequency.  For  99%  of  the  data,  asc  is  less  than  4.66  Hz,  1.38  Hz,  0.22  Hz,  respectively,  at 
VHF,  UHF,  L-band.  Clearly,  this  spectral  spread  decreases  to  very  small  values  as  the  frequency 
increases.  To  get  a  worst-case  (largest)  value  of  asc >  substitute  the  smallest  value  of  ro  (0.024  s) 
into  Eq.(38)  to  obtain  9.38  Hz  for  the  three  frequencies.  These  values  are  then  used  in  Eqs. (31) 
and  (34). 
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Fig.  7  —  Cumulative  distribution  of  the  clutter  spread  Osc  caused  by  ionospheric 
scintillation  for  the  1977  Ancon  satellite  passes  at  VHF,  UHF,  and  L-band 
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4.3  Computation  of  a  Mainbeam’s  Clutter-Doppler  Spectrum 

To  determine  the  total,  clutter- Doppler  spread  of  a  mainbeam,  consider  the  spectra  Scl  and 
Scl,a ,  obtained  by  substituting  Eqs.(31)  and  (34)  into  Eq.(25).  These  spectra  are  graphed  for 
nominal  values  of  rc,  h,  fr,  po  (1  ps,  1030  km,  577  Hz,  1°)  and  for  five  sets  of  (vr,4>a,6 0)  at  the 
three  frequencies.  Recall  that  rc,  h,  fr,  and  po  are  the  compressed  pulsewidth,  satellite’s  altitude, 
PRF,  and  one-half  beamwidth.  Three  of  the  sets  are  taken  from  the  Ancon  data  (the  dark  boxes 
in  Table  4  and  lines  1  to  3  of  Table  5).  The  remaining  two  sets  are  (0,21.731°, 282.61°)  and 
(7789.230 m/s, 58.297°, 90°).  Approximate  values  of  vr  for  the  data  (Sets  1  to  3)  are  generated 
from  Eq.(35),  with  an  exception  discussed  below,  while  Sets  4  and  5  are  chosen  so  that  &pm  is 
extremized  for  a  1030-km  orbit. 


Table  5  —  Five  Sets  of  Data 


Since 

v\  =  Yr  •  Yr  =  v2p  +  v\  -  2(vp  •  ve), 

the  maximum  value  of  vr  is  attained  when  ve  and  vp  are  antiparallel.  Hence 

VT  <  Vp  +  Ve. 

By  Eqs.  (37)  and  (39),  an  upper  bound  for  Opm  is 


7T 

720 


I  8  2  ft ,  x 
hT20T^  +  ^’ 


(39) 

(40) 


where  the  speeds  ve  and  vp  are  calculated  in  Appendix  A. 

Recall  from  Eq.(  15)  that  only  angles  (f>  out  to  the  optical  horizon  59.425°  are  considered;  that  is, 
<j>a  <  58.425°  for  a  2°  antenna  beamwidth  (2po).  Consequently,  for  a  1030-km  orbit,  |  sind>„|  <  0.852 
and 

km  km 

vr  <  0.852(up  +  »e)  ^  (0. 852 )(7. 78923—  =  6.6364—.  (41) 

Therefore,  6.6364  km/s  is  the  maximum  achievable  value  of  vr  for  this  orbit;  however,  observe  that 
this  value  is  exceeded  in  Sets  1  to  3  of  Table  5.  This  merely  indicates  that  the  approximation 
(Eq.(37)),  used  to  calculate  vT  for  the  data,  is  an  overestimate.  In  lieu  of  no  velocity  measurement 
and  incomplete  orbital  data,  the  inexactitude  (upper  bound)  in  the  calculation  of  vr  is  accepted. 
The  exception  is  Set  2,  where  the  value  9.5471  km/s  of  Table  4  is  replaced  with  vp  +  ve  =  7.7892 
km/s  to  compensate  for  the  poor  approximation  of  vr  •  Rs  by  vp  •  R5  in  this  case  (see  Appendix 
A).  Correspondingly,  the  values  of  a  pm  of  Table  4  are  changed  to  18.43  Hz,  55.30  Hz,  and 
165.88  Hz  at  VHF,  UHF,  L-band,  respectively.  Note  that  these  modified  values  are  still  the  maxi¬ 
mum  <7 pm  achieved  over  the  data. 
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The  selection  of  577  Hz  and  1  fis  for  the  PRF  and  compressed  pulsewidth  correspond  to  1733 
(Nr)  range  cells  per  unambiguous  range  interval  (Eq.(14)).  The  number  of  unambiguous  range 
intervals  within  the  mainbeam  are  given  in  Table  5.  Set  1  corresponds  to  the  data  set  for  which 
the  Wideband  boresight  pointed  nearest  the  optical  horizon  ( <j>  ~  59.425°).  Hence  the  number  of 
unambiguous  range  intervals  for  set  1  is  the  largest  of  all  the  Ancon  data.  Sets  2  and  3  correspond 
to  the  maximum  and  minimum  values  of  opm  for  the  data  (Table  4).  The  last  two  sets  are  the 
minimum  (0  Hz,0  Hz,0  Hz)  and  the  maximum  (36.98  Hz,  110.93  Hz,  332.79  Hz)  values  that  a  pm 
can  attain  for  VHF,  UHF,  and  L-band  and  for  an  arbitrary,  1030-km  orbit. 
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f  (Hz) 

Fig.  8  —  Exact  and  approximate  mainbeam  clutter-Doppler  spectra  for  set  #1  at  VHF,  UHF,  and  L-band. 
Both  spectra  are  normalized  to  the  maximum  value  of  the  exact  spectrum.  Twice  the  horizontal  distance 
between  the  frequency  of  the  maximum  value  of  each  spectrum  and  the  vertical  lines  to  both  curves 
correspond  to  the  clutter  spread  for  each  function. 


21 


Normalized  Mainbeam  Normalized  Mainbeam  Normalized  Mainbeam 

Clutter-Doppler  Spectrum  Clutter-Doppler  Spectrum  Clutter-Doppler  Spectrum 


MOKOLE 


Fig.  9  -  Exact  and  approximate  mainbeam  clutter-Doppler  spectra  for  set  #2  at  VHF,  UHF,  and  L-band. 
Both  spectra  are  normalized  to  the  maximum  value  of  the  exact  spectrum.  Twice  the  horizontal  distance 
between  the  frequency  of  the  maximum  value  of  each  spectrum  and  the  vertical  lines  to  both  curves 
correspond  to  the  clutter  spread  for  each  function. 
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Fig  10  -  Exact  and  approximate  mainbeam  clutter-Doppler  spectra  for  set  #3  at  VHF,  UHF,  and  L-band. 
Both  spectra  are  normalized  to  the  maximum  value  of  the  exact  spectrum.  Twice  the  horizontal  distance 
between  the  frequency  of  the  maximum  value  of  each  spectrum  and  the  vertical  lines  to  both  curves 
correspond  to  the  clutter  spread  for  each  function. 
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Fig.  1 1  —  Exact  and  approximate  mainbeam  clutter-Doppler  spectra  for  set  #4  at  VHF,  UHF,  and  L-band. 
Both  spectra  are  normalized  to  the  maximum  value  of  the  exact  spectrum.  Twice  the  horizontal  distance 
between  the  frequency  of  the  maximum  value  of  each  spectrum  and  the  vertical  lines  to  both  curves 
correspond  to  the  clutter  spread  for  each  function. 
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Fig.  12  -  Exact  and  approximate  mainbeam  clutter-Doppler  spectra  for  set  #5  at  VHF,  UHF,  and  L-band. 
Both  spectra  are  normalized  to  the  maximum  value  of  the  exact  spectrum.  Twice  the  horizontal  distance 
between  the  frequency  of  the  maximum  value  of  each  spectrum  and  the  vertical  lines  to  both  curves 
correspond  to  the  clutter  spread  for  each  function. 
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Table  6  —  Clutter- Doppler  Spectral  Spreads  and  Maxima 


Set 

Freq 

OpM 

(Hz) 

OCL 

(Hz) 

OCL 

(Hz) 

OCL, A 

(Hz) 

arg{max(5ci)} 

(Hz) 

max(ScL) 

(xlO-18) 

1 

VHF 

9.95 

13.71 

22.39 

28.02 

-2.123 

1.77651 

1 

UHF 

29.85 

31.30 

50.25 

79.83 

-6.576 

0.62352 

1 

L-band 

89.54 

90.03 

142.96 

238.00 

-19.716 

0.20909 

2 

VHF 

18.43 

20.70 

21.25 

22.87 

0.104 

6.41250 

UHF 

55.29 

56.09 

57.63 

63.22 

0.316 

2.31994 

L-band 

165.88 

166.15 

171.03 

187.87 

0.952 

0.78098 

VHF 

4.55 

10.47 

22.85 

27.51 

1.230 

2.34508 

UHF 

13.65 

16.59 

29.29 

78.25 

4.010 

0.82969 

L-band 

40.94 

42.01 

61.49 

233.60 

12.306 

0.27828 

4 

VHF 

0.00 

9.43 

9.48 

9.48 

0.000 

15.06030 

4 

UHF 

0.00 

9.43 

9.48 

9.48 

0.000 

15.06030 

4 

L-band 

0.00 

9.43 

9.48 

9.48 

0.000 

15.06030 

5 

VHF 

36.98 

38.16 

36.42 

37.93 

0.016 

1.31096 

5 

UHF 

110.93 

111.33 

106.21 

110.81 

0.000 

0.44956 

5 

L-band 

332.79 

332.92 

317.53 

331.54 

0.083 

0.15034 

A  natural  question  is  whether  the  spectra  Scl  and  Scl, A  are  Gaussian.  These  spectra  are 
graphed  in  Figs.  8  to  12.  Let  / m  and  / ma  be  the  maxima  associated  with  Scl  and  Scl, a, 
respectively.  Although  Gaussian  fits  to  the  data  are  not  shown,  sets  2  and  5  are  nearly  Gaussian, 
sot  4  is  Gaussian,  and  sets  1  and  3  are  not.  Clearly,  set  4  is  Gaussian  because  vr  =  0;  hence  the 
exponential  function  is  identical  for  every  summand  of  Eq.(25).  When  normalized,  a  finite  sum 
of  Gaussians  of  equal  standard  deviations  and  zero  means  is  a  zero-mean  Gaussian  with  the  same 
standard  deviation,  evidenced  by  the  fact  that  ocl  —  &cl  in  Table  6. 

Sets  1  and  3  are  also  not  symmetric  about  f\j,  since  more  area  beneath  Scl  lies  to  the  right 
of  f\i  than  to  the  left  of  it  (Figs.  8  and  10).  Nonetheless,  the  standard  deviations  &cl  and  ocl,a 
are  defined  to  be  the  frequencies  for  which  the  areas  under  Scl  and  Scl, a  within  the  intervals 
[/a/,  f\l  +  °cl\  and  [ /mAi/ma  +  °cl,a]  are  0.6826  of  the  area  to  the  right  of  / m  and  /ma-  As 
the  spectra  become  more  symmetric  about  the  maxima,  this  definition  agrees  more  closely  with 
the  distance  from  the  1-sigma  points  to  the  center  (maximum)  of  the  fitted  Gaussian  probability 
density  function.  The  graphs  of  the  exact  spectra  are  close  enough  to  being  symmetric  to  yield 
reasonable  values  for  ocl-  Intuitively,  one  expects  the  summation  of  many  Gaussians  with  different 
means  to  result  in  spectral  spreading  (ocl  >  ocl)-  This  is  true  except  for  set  5,  where  &cl 
is  underestimated  because  Scl  matches  a  Gaussian  until  j/|  =  700  Hz,  while  the  exact  spectrum 
decreases  more  rapidly  to  zero  than  the  fitted  Gaussian  for  |/|  >  700  Hz.  In  addition,  the  difference 
between  ctcL  and  ocl  for  vT  ^  0  is  more  pronounced  as  the  frequency  increases. 

Observe  that  the  maximum  /m,  which  is  taken  as  the  center  of  the  spectrum  Scl,  is  not 
approximately  zero  in  all  cases.  In  particular,  the  maximum  in  sets  1  and  3  can  be  several  Hz  (as 
much  as  2,  6,  and  19' for  VHF,  UHF,  and  L-band,  respectively)  away  from  0  Hz.  This  implies  that 
the  TACCAR-like  correction  (Eq.(7))  to  the  Doppler  shift  does  not  necessarily  center  the  spectrum 
at  0  Hz.  Table  6  also  indicates  that  the  maximum  value  of  the  spectrum  decreases  with  increasing 
frequency  at  a  given  nonzero  vr. 

Upon  analyzing  columns  5  and  6  of  Table  6,  it  is  clear  that  (Jcl,a  is,  at  best,  a  reasonable 
approximation  of  ocl  (Sets  2  and  5)  for  nonzero  vr.  In  fact,  the  difference,  ocl,a  ~  ®cl,  is 
nonnegative  in  all  instances,  increases  with  frequency,  is  substantial  for  sets  1  and  3,  and  is  zero 
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for  set  4  (since  B  =  0).  Therefore,  unless  vr  is  small  enough,  the  representation  of  the  spectrum 
given  by  Eqs.(31)  and  (32)  should  be  used. 

Recall  that  apj(v)  is  1  m/s  and  a  sc  is  9.38  Hz.  Since  a(f)  =  2fta(v)/c  relates  the  frequency- 
and  velocity- dependent  spreads,  api(f)  is  0.92  Hz,  2.76  Hz,  and  8.27  Hz  at  VHF,  UHF,  and  L-band, 
respectively.  Depending  on  the  antenna’s  position  in  a  1030-km  orbit,  apM  ranges  from  0  Hz  to 
369  Hz.  Clearly,  apj-  <  a  so  for  these  calculations.  At  VHF,  as c  >  &in  and  °in  contributes  very 
little  to  a  cl  in  the  absence  of  apM-  However,  as  the  frequency  increases  to  L-band,  ap;  and  asc 
are  nearly  the  same  and  contribute  almost  equally  to  acL  when  apM  ~  0.  On  the  other  hand,  a  pm 
can  be  large  compared  to  asc  and  &IN,  in  which  case  a  pm  —  a  cl • 

As  indicated  previously,  one  expects  spectral  spreading  from  integration  over  the  antenna’s 
mainbeam;  that  is,  a  cl  is  larger  than  ocl-  The  spreading  is  a  minimum  (act  ~  &CL  =  0.0538  Hz) 
when  apM  is  zero  and  increases  with  frequency,  except  for  the  noted  aberration  of  set  5. 

5.0  SUMMARY 

Two  expressions  for  the  clutter-Doppler  spectrum  within  a  range  cell  and  for  the  total  spectrum 
over  a  narrow,  Gaussian  mainlobe  are  derived  for  an  SBR.  These  expressions  are  applicable  to  an 
SBR  in  a  1030-km  orbit  with  a  carrier  frequency  between  100  MHz  and  1300  MHz  during  severe 
ionospheric  conditions.  In  all  situations,  it  is  more  desirable  to  use  the  more  exact  of  the  two 
expressions,  except  when  the  speed  of  the  satellite  relative  to  the  speed  of  the  surface  clutter  patch 
is  “near”  zero.  The  depth  of  this  study  is  insufficient  to  determine  a  quantitative  definition  of 
nearness. 

A  total  spectrum  is  not  necessarily  symmetric  about  its  maximum;  however,  each  spectrum 
is  sufficiently  symmetric  to  define  a  center  (maximum)  and  a  spread  about  the  maximum  that  are 
analogous  to  the  mean  and  standard  deviation  of  a  Gaussian  probability  density  function.  At  VHF, 
UIIF,  and  L-band,  the  center  of  the  mainbeam ’s  spectrum  can  be  offset  from  zero  by  as  much  as 
2.12,  6.58,  and  19.72  Hz,  respectively.  Since  these  calculations  include  a  TACCAR-like  correction 
to  the  Doppler  shift,  such  a  correction  does  not  compensate  completely  for  the  translation  of  the 
spectrum  caused  by  platform  motion.  If  the  offsets  from  0  Hz  are  significant  to  an  SBR  system, 
then  an  additional  correction  is  required  to  center  the  spectrum  at  zero;  however,  making  this 
correction  is  limited  by  the  ability  to  perform  such  signal  processing. 

Calculations  of  asci  vcl,  and  a  cl  are  made  for  the  minimum  value  of  tq  (24  ms)  over  the 
Wideband  data  collected  at  Ancon  (Peru)  from  22  satellite  passes  in  1977.  This  minimum  value 
is  typical  of  severe  scintillation  conditions  and  yields  a  a  sc  of  9-38  Hz.  For  the  data,  the  intrinsic 
clutter-Doppler  spread  a  in  ranges  from  0.92  Hz  to  8.27  Hz.  Depending  on  the  orbital  position, 
a  pm  takes  on  values  between  0  Hz  and  369  Hz.  When  apM  is  near  zero  and  a  in  is  small  compared 
to  asc ,  for  example  at  VHF,  acL  —  &9C-  Situations  also  exist  where  asc  and  cr^v  are  comparable 
for  small  a  pm  •  On  the  other  hand,  as  a  pm  increases,  it  dominates  a  sc  and  a  in  so  that  a  cl  —  &PM- 
All  of  these  situations  occur  or  can  occur  for  a  satellite  in  some  1030-km  orbit.  Thus,  in  some  cases, 
eliminating  the  clutter-Doppler  spread  caused  by  platform  motion  with  techniques  like  DPCA  result 
in  a  negligible  decrease  in  the  spectral  spreads,  whereas  in  other  instances,  the  decrease  is  dramatic. 
At  best,  one  can  reduce  acL  and  acL  only  to  y/ajN  +  a2^  ~  9.43  Hz  and  9.48  Hz,  respectively,  in 
severe  scintillation.  The  preceding  two  limits  provide  lower  bounds  on  the  achievable  reduction  in 
clutter-Doppler  spread.  To  obtain  a  more  quantitative  idea  of  the  ramifications  of  these  bounds, 
one  must  implement  a  system  architecture  and  analyze  improvement  factors,  clutter  cancelers,  and 
the  like. 
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Appendix  A 

APPROXIMATION  OF  THE  DOPPLER  SHIFT 


The  Doppler  shift  of  a  scatterer  at  point  Ps  in  the  direction  Rs  is 

f D  =  ~{'iT  •  Rs), 

where  vr  =  vp  -  ve.  Two  approximations  can  be  made  in  Eq.(4):  replacement  of  vr  •  Rs  with 
vp  •  Rs  or  7r  with  7 s-  This  appendix  addresses  the  accuracy  of  these  replacements  for  a  circular 
orbit  at  1030  km. 

The  speed  of  the  rotating  point  Ps  is 


Ve 


|ve|  =  |w  x  raj  =  |(weK)  x  ij 

We(— +  =  W ey/r^x  +  r]Y  =  ueRe 


cos  L 


(Al) 


where 


Jr  rad 
Ue  ~  (12)(3600)  s  ’ 

fgX  —  Re  COS  L  COS  £, 

rgY  =  RecosLsin£.  (A2) 


The  angles  L  and  f  are  the  latitude  and  longitude,  respectively,  of  the  ground  station  (see  Fig. 
1(c)),  and  Re  (6378.145  km)  is  the  radius  of  the  Earth.  Therefore,  ve  =  ueRe  cos  L,  which  is 
independent  of  the  longitude  of  the  ground  station.  For  the  site  at  Ancon,  Peru,  the  latitude  is 
-11.767°.  Hence  ve  ~  454.0843718  m/s. 

According  to  (Eisele  and  Nichols  1976,  p.12),  the  orbital  speed  of  a  satellite  in  m/s  is 


Vp  ~  (4.2685)(1852) 


■R 


Rj-r  ~  7335.145367 -, 
+  h  s 


(A3) 


where  h  (1030  km)  is  the  satellite’s  altitude. 

Comparing  the  two  speeds  yields  ve  ~  0.0619053  t/p;  that  is,  ve  is  about  6.19%  of  vp  for  a 
1030-km  orbit.  Thus  for  cos  75  /  0, 


¥r  •  Rs  =  v  •  Rs  -  ve  •  Rs 


^  v. 


cos  7s  -  0.0619053  cos  if>& 


COS  tA.c 

0.0619053 - — 

cos  7s 


(A4) 


where  75  and  V^eS  are  the  angles  between  vp  and  Rs  and  ve  and  Rs,  respectively.  The  constant  D 
is  the  second  bracketed  expression  in  the  next  to  last  line  of  Eq.(A4).  Whether  the  term 


cos  d>„<f 

0.0619053 - — 

cos  7s 


(AS) 
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can  be  neglected  depends  on  the  satellite’s  orbit  and  the  direction  of  antenna’s  boresight.  When 
the  magnitude  of  this  expression  is  large  (>  0.1),  vr  •  Rs  cannot  be  accurately  approximated  by 
vp  •  Rs-  For  example,  if  7 s  is  near  jt/2  and  i^es  is  not,  then  this  quantity  gets  very  large.  In  fact,  it 
approaches  infinity.  If  9s  =  x/2,3x/2  (near  broadside)  or  <£s  =  0  (nadir),  then  vr  •  Rs  =  -ve  -  Rs- 
Therefore,  the  component  of  fp  from  the  platform  speed  is  zero,  when  the  mainbeam  points 
broadside  or  in  the  nadir  direction. 

At  the  equator,  ve  is  perpendicular  to  the  nadir  direction  and  points  along  the  equator.  Anal¬ 
ysis  of  the  geometry  leads  to 


*  .  Re 

—  -  arcsin  — - - 

2  \  Re  +  h 


<**<£  + arcsin 


or 


Re 

Re  4"  h. 


<  COS  IpeS  < 


Re 

Re  +  h 


(A7) 


for  radars  that  only  look  out  to  the  optical  horizon. 

To  determine  the  values  of  D  for  the  Wideband  data,  ipes  must  be  expressed  in  terms  of 
measurable  quantities.  The  radar  site  measures  the  position  p  of  the  satellite  relative  to  the  site’s 
position  in  topocentric  coordinates.  In  Fig.  1,  p  corresponds  to  —  R5.  In  particular, 


P  =  ps  S  +  Pe  E  +  pz  Z, 


(A8) 


where 


ps  =  - p  cos  El  cos  Az, 

Pe  —  P  cos  El  sin  Az, 
pz  =  psin  El. 

(A9) 

Therefore, 

Rs  =  (cos  El  cos  Az)S  +  (—  cos  El  sin  Az) E  +  (-  sin  El) Z 

(A10) 

and 

ve  •  Rs  =  (veE)  •  Rs  =  Ve  COS  TpeS- 

(All) 

Combining  Eqs.  (A10)  and  (All)  yields 

cos  ipes  =  -  cos  El  sin  Az. 

(A12) 

Therefore, 

„  cos  El  sin  Az 

D  =  1  +  0.0619053 - r — 

cos  9$  sin  <ps 

(A13) 

To  use  the  Wideband  data,  set  <f>s  =  4>a  and  9s  =  9a.  Since  El ,  Az,  <f>a,  and  9a  are  measured  and 
tabulated  in  Table  2,  the  values  of  D  and  the  absolute  error,  |(Z)  -  \)/ D\,  are  computed  in  Table 
Al.  The  error  is  good  only  for  the  last  three  rows  of  Table  A1  and  is  less  than  10%  for  all  but  one 
the  data  sets,  where  it  is  a  substantial  23.77%.  Unfortunately,  this  is  one  of  the  sets  analyzed  in 
Table  5.  In  fact,  the  errors  are  2.87%,  3.89%,  and  23.77%  for  Sets  1,  2,  and  3,  respectively,  of  Table 
5. 


NRL  REPORT  9327 


Table  Al  —  Values  of  D  and  the  Absolute  Error 
for  an  Orbit  Relative  to  Ancon 


Az 

(Deg) 

El 

(Deg) 

<t>a 

(Deg) 

0a 

(Deg) 

D 

Absolute  Error 
(%) 

23.6 

10.4 

58.297 

346.34 

1.0295 

2.87 

25.5 

15.8 

56.503 

344.68 

1.0319 

3.09 

28.1 

22.4 

53.516 

342.25 

1.0352 

3.40 

32.3 

30.9 

48.752 

338.48 

1.0405 

3.89 

39.8 

42.0 

41.427 

331.68 

1.0506 

4.82 

56.1 

55.7 

31.188 

317.52 

1.0758 

7.05 

96.6 

65.8 

21.731 

282.61 

1.3119 

23.77 

124.1 

64.7 

24.411 

232.50 

0.9129 

9.54 

157.9 

52.8 

35.007 

207.69 

0.9723 

2.85 

171.3 

39.7 

44.066 

197.24 

0.9892 

1.09 

177.8 

29.3 

50.231 

191.86 

0.9972 

0.28 

-178.5 

21.4 

54.182 

188.64 

1.0019 

0.19 

-176.2 

15.1 

56.632 

186.50 

1.0048 

0.48 

To  avoid  approximations  to  the  Doppler  shift,  exact  expressions  must  be  employed  when 
analyzing  vr  •  Rs.  One  major  approximation,  which  greatly  simplifies  the  analysis,  is  substituting 
7s  for  7r.  To  obtain  7r  exactly,  express  Rs  and  yr  in  the  vehicle- centered  frame  ( x,y,z ): 

Rs  =  *SX  +  ysy  +  zSZ  and  vr  =  vrxx  +  vryy  +  vrzz, 


where 

xs  =  Rs  cos  0S  sin  <f>s 
Vs  =  RssmOs  sin  <j>s  and 
zs  =  Rs  cos  4>s 
Equating  both  sides  of 


vrx  —  vr  cos  0r  sin  <f>r 
vTy  =  vr  sin  0r  sin  <f>r 

Vrz  =  Vr  COS  4>r- 


Vr  COS7r  —  Yr  '  R"S  —  D  (^n %S  "b  ^ ryVS  "1"  ^rz^s) 

Rs 


leads  to 


cos  7 r  =  cos  0r  sin  4>t  cos  0S  sin  4>s  +  sin  0r  sin  <f>r  sin  Os  sin  <f>s  +  cos  <t>r  cos  <f>s.  (A14) 

The  computation  of  <f>r  and  0T  requires  precise  information  about  both  the  position  and  the  velocity 
of  the  satellite  with  respect  to  the  three  coordinate  frames  that  are  mentioned  in  Sections  1  and  2. 
Lastly,  one  more  expression  for  vr  •  Rs  is  derived.  By  Eqs.(A8)  and  (A9), 

P  =  Ps§  +  PeE  -I-  pzZ,  (A15) 

where 

ps  =  -p  cos  El  cos  Az  -f  pEl  sin  El  cos  Az  +  pAz  cos  El  sin  Az , 

Pe  =  P  cos  El  sin  Az  -  pEl  sin  El  sin  Az  +  pAz  cos  El  cos  Az, 

pz  =  p  sin  El  -(-  pEl  cos  El.  (A16) 

Since  yr  =  p, 

Yr  •  R s  =  p-  =  +  PeV  +  PzZ)  =  -(>  =  Rs,  (A17) 

after  some  algebra.  When  applying  the  Wideband  data,  Rs  =  Ra. 
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Appendix  B 

APPROXIMATION  OF  AVERAGE  DOPPLER  OVER  NARROW  MAINBEAMS 


An  approximate  expression  for  the  average  Doppler  over  a  narrow  mainbeam  is  derived  for 
isotropic,  homogeneous  clutter.  By  Eq.(5), 

j  ,,  ff  fD(Ro,<t>,0)R%sin<l>d4>d6 

/d=  A^(AjjJ  - ff  Rl  sin  J  d4M - •  (B1) 


Assuming  that  vT  is  essentially  constant  over  fl  and  substituting  /d(#o,<£,  0)  =  (2  ftvr/c)  sin  cf>  cos  6 
into  Eq.(Bl)  yields 

//  cos  6  s\tl2 4>  d(j>dO 

jo  ~  2ftVr  n  .  ,  (B2) 

c  J  J  sin  4>  d(f>d6 

n 

where  4>  and  8  run  over  all  scatterers  in  the  mainbeam  and  the  subscript  S  of  Eq.(4)  is  dropped. 

The  geometry  of  Q  clearly  determines  the  evaluation  of  Eq.(B2).  However,  for  narrow  main- 
beams,  the  right  side  of  Eq.(B2)  reduces  to  the  same  approximate  formula  for  small  (c  and  6  are 
small)  rectangular  (fli)  and  elliptical  (O2)  neighborhoods  about  (<j>a,0a): 

fil  =  |(<M)  :  4>a  -  f  <  <t>  <  4>a  +  C  &  0a  -  0  <  0  <  0a  + 

=  <B3) 

Since  evaluation  of  the  integrals  over  fii  in  Eq.(B2)  is  straightforward,  the  approximation  is  derived 
for  fl2-  VVith  the  change  of  variables  from  to  (u,v), 

(f>  =  tu  COS  V  +  <f)a 

(B4) 

SI2  and  the  areal  element  d<f)d9  become  {(u,u)  :  0  <  u  <  1  &  0  <  v  <  27r}  and  (Sududv, 

respectively. 

First,  consider  the  integral  in  the  numerator  of  Eq.(B2). 

Jj  cos  0sin2  4>d<t>dd  =  —  J  J  ujcos(u£  sin  v  +  9a)  |l  -  cos(2uecosv  +  2<f>a)  |  dudv.  (B5) 


In  the  two  terms  of  the  inner  integral,  use  trigonometric  identities,  integrate  by  parts,  and  apply 
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the  small-angle  approximations,  sin  z  ~  z  and  cos  z  ~  1,  to  obtain 
J  ucos(u6  sin  v  +  9a)  cos(2uccos  v  +  2<f>a)  du  =  ^  J  u  jcos[ti(£ sin  v  -  2c cos  v)  +  ( 8a  -  2 </>„)] 

+  cos[u(£sint>  +  2c  cos  v)  +  (8a  +  2<£a)]  |  du 

if  /sin [u(£ sin  v  -  2c  cos  v)  +  ( 6a  -  2<j>a)] 

2  1  V  S  sin  v  —  2c  cos  v 

sin  [u(6  sin  v  +  2c  cos  v)  +  ( 8a  -(-  2<£a)]  \ 


1  J  (£sinu 

H 


6  sin  v  +  2e  cos  v 

/ 

cos  [u(£ sin  v  —  2c  cos  v)  +  ( 8a  -  2<f>a)] 
(8  sin  v  -  2e  cos  t>)2 

cos  [u  (£  sin  u  2c  cos  u)  +  ( 8a  +  2  <j>a )  ] 
(8  sin  v  +  2c  cos  v )2 

-  2c  cos  v ) cos (9 a  -  2<f>a)  +  sin(0a  -  2 <f>a) 

6  sin  v  —  2c  cos  v 


+ 


(6  sin  v  +  2c  cos  n)  cos(ga  +  2 <j>a)  -|-  sin(fla  -f  2 <f>a) 
8  sin  v  -)-  2c  cos  v 

sin(0a-2<£a)  sin(0a  +  2<£a)  1 


6  sin  v  -  2c  cos  v  8  sin  v  +  2c  cos  v  j 

=  cos  8a  cos(2<£0). 

Similar  manipulations  lead  to 

I  ucos(u8  sinv  +  9a)du  ~  cos0a. 

Jo 

Substituting  Eqs.(B6)  and  (B7)  into  Eq.(B5)  yields 

JJ  cos  9  sin2  4>dcf>d8  ~  j  J  cos9a  [l  -  cos(2</>a)]  dv 
n2  0 

=  27rc 6  cos  8 a  sin2  <f>a. 

Evaluation  of  the  denominator  of  Eq.(B2)  is  much  simpler;  that  is, 

j  J  sin  <j>  d<t>d9  =  (S  j  u  sin(cu  cos  v  +  4>a )  dudv 
n2  0  0 

-  8  [2* l  «cos(cucosv  +  <j>a)  sin(cu cos v  +  4>a)' 
Jo  \  c  cos  V  (c  cos  v)2 

r2ir 

~  (6  I  sin  <f>a  dv  =  27rc6  sin  <f>a. 

Jo 

Combining  Eqs.(B2),  (B8),  (B9)  yields  the  desired  result, 

jd  - - COS  8a  sin  <pa. 


dv 


(B6) 

(B7) 


(B8) 


(B9) 

(BIO) 
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RESIDUAL  DOPPLER 


By  Eqs.(7)  and  (23), 

F'n^mnj)  =  /o(^mn  +  <j>a,9  +  0a)  =  {  cos  (0  +  6a)  sin(4>mn  +  4>a)  -  cos  6  a  sin  <j>a  | 

cos  9  sin  <f>mn  cos  6a  cos  <f>a  -f  cos  0  cos  <^mn  cos  0a  sin  <f>a 
—  sin0sin<£mflsin0acos0a  -  sin  0  cos  <}>mn  sin  9a  sin  (j>a  —  cos0asin</>a  j.  (Cl) 

Suppose  the  mainbeam  is  narrow;  that  is,  0  and  <t>mn  are  less  than  or  equal  to  2°.  A  first-order 
approximation  to  F'D  is  obtained  by  keeping  only  first-degree  terms  of  9  and  <j>mn  in  the  products 
of  the  Maclaurin  expansions  of  the  trigonometric  functions  that  depend  on  them.  In  particular, 

cos0sin^mn  ~  4>mn, 
cos0cos4>mn  ~  1, 
sin0sin<^mn  ~  0, 

sin^cos^mn  ~  9.  (C2) 

Hence 

F'D(4>mn,9^  cos  0  a  COS  <l>a  -  0  sin  9  a  sin  <f)ay  (C3) 
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Appendix  D 

INTEGRALS  INVOLVING  THE  ANTENNA  POWER  GAIN 


Two  points  are  discussed  in  this  appendix.  First,  the  relationship  7?  =  poy/%!  ln(20)  is  estab¬ 
lished.  Second,  the  accuracy  of  replacing  Eq.(31)  by  Eq.(34)  is  analyzed. 

The  mainlobe  of  the  Gaussian  gain  of  Eq.(31)  is  identified  with  the  angle  p0,  where  2po  is  the 
mainlobe’s  beamwidth  (Fig.  Dl).  The  selection  of  po  depends  on  the  desired  beamwidth.  This 
report  assumes  that  95%  of  the  power  lies  within  the  region  Z  -  {(<M)  :  0  <  <j>2  +  62  <  pi}.  This 
implies  that 


Fig.  Dl  —  The  mainlobe  of  the  two-way 
gain  of  a  Gaussian  antenna  pattern 


With  the  change  of  variable  4>  =  pcosS  and  6  =  psintf,  the  preceding  equation  becomes 

r2ir  r  00 

/  9oe~*P* ^ pdpdti  =  0.95  /  Poe_8p  ^  pdpdS, 

Jo  Jo 
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or  equivalently, 


V  =  Po\ 


ln(20) 


=  PoM, 


where  /i  ~  0.611936708. 

Next,  consider  the  integral  of  Eq.(29), 


exp(  2^+^)  ) M-  (D2) 

Since  the  integrand  is  bounded  above  by  \G\A,  the  integral  in  Eq.(D2)  has  an  upper  bound  that  is 
the  integral  of  \Q\ 4  over  the  same  range  of  integration.  Therefore,  if 


r°° 

/  /  I G(*mn,S)\4d9, 

J  _  J  —  OO 


_  4  (  [f  -  t 

L_  “p( — 2(? 


(f  - 


tin  +  asc) 


/OO 

-OO 


X  exp  I  - 


(f  - 

2Kv  +  °lc) 


By  Eq.(30)  and  dlnm  =  -1?^, 


=  -l|,  (D4) 

where  the  Normal  Probability  Function  P  is  (Abramowitz  and  Stegun  1968) 

P{*)  =  I  -4=e_‘2/2  dt-  (D5) 

J  —  00  V 

The  integral  of  Eq.(D4)  is  evaluated  over  the  vertical  slice  depicted  in  Fig.  D2.  By  Eq.(D4)  and 
l>mn  =  M  ~  4>lmi 

J*t  \G{<f>mn,$)\*  dd  =  ?02^e-w™/r,J|2P^^  -  -  l|. 

To  determine  the  accuracy  of  Eq.(D3),  analyze  the  ratio 


W(4>m „)  = 


=  T-,SV',-ft 
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As  the  ratio  decreases,  the  accuracy  of  Eq.(D3)  improves.  Clearly,  Eq.(D6)  depends  on  how  close 
4>mn  is  to  po-  In  particular,  by  applying  Table  26.1  of  (Abramowitz  and  Stegun),  W  attains  its 
maximum  (minimum)  value  of  unity  (0.0144)  when  4>  =  ±Po  (0).  The  actual  integral  is  98.56%  of 
the  approximation  for  <j>  =  0,  which  is  reasonably  good;  however,  the  approximation  is  the  worst  it 
can  be  at  <f>  =  ±po . 


Fig.  D2  —  The  region  that  is  subtended  by  the  Gaussian 
pattern’s  roainlobe  in  the  coordinates  ^  and  6 


Let  <j>  —  Apo,  —  1  <  A  <  1,  correspond  to  an  arbitrary  <£-slice  in  the  mainbeam  (Fig.  D2).  Then 
the  argument  of  P  becomes  4\/l  —  A2/p.  Since  W(4>)  =  0.0930  for  A  =  0.7273,  the  approximation  is 
still  fair  for  all  \<f>rrm\  <  0.7273po-  Although  the  approximation  is  not  good  for  \<j>rrm\  6  [0.7273po)Po]> 
the  amplitudes 

(D7) 

of  these  contributions  to  the  clutter-Doppler  spectrum  are  <  0.20502  of  the  contribution  at  <f>  =  0 
(see  Fig.  D3  and  Table  D1  for  evaluations  of  exp(-8<^„/?72)).  Consequently,  they  won’t  contribute 
as  heavily  to  Scl{})-  In  particular,  the  spectral  contributions  for  strips  near  the  heel  and  toe  of  the 
mainbeam  may  be  ignored.  To  determine  whether  the  contributions  in  the  transitional  region,  from 
0.7273po  to  a  value  a  bit  less  than  po>  one  must  analyze  them  for  each  power  gain  on  a  case-by-case 
basis. 
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Fig.  D3  —  The  exponential  term  of  the  amplitude  in 
Eq.(D7)  vs  the  fraction  A  of  the  half-beamwidth  p0 


Table  D1  —  Exponential  vs  A 


A 

exp(-  S<))2/t]2) 

0.0000 

1.00000 

0.2500 

0.82925 

0.5000 

0.47287 

0.7273 

0.20502 

0.7500 

0.18543 

0.8000 

0.14701 

0.8500 

0.11482 

0.9000 

0.08834 

0.9500 

0.06696 

1.0000 

0.05000 

REFERENCE 

M.  Abramowitz  and  I.  A.  Stegun,  ed.,  Handbook  of  Mathematical  Functions  (Dover  Publications, 
Inc.,  New  York,  1968),  p.  970. 


39 


